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ABSTRACT 

Context. The current detectors of gamma-ray emission specially at TeV energies have too poor resolution to determine 
whether this emission is produced in the jet or in the core, particularly of low luminous, non-blazar AGNs (like radio 
galaxies). In recent works it has been found that the power released by events of turbulent fast magnetic reconnection 
in the core region of these sources is more than sufficient to reproduce the observed gamma-ray luminosities. Besides, 
3D MHD simulations with test particles have demonstrated that a first-order Fermi process within reconnection sites 
with embedded turbulence results very efficient particle acceleration rates. 

Aims. We computed here the spectral energy distribution (SED) from radio to gamma-rays of the radio galaxies for 
which energy emission up to TeVs has been detected (namely, M87, Gen A, Per A, and IG 310). 

Methods. For this aim, we employed the acceleration model above and considered all the relevant leptonic and hadronic 
loss processes around the core region of the sources. 

Results. We found that the calculated SEDs match very well specially with the VHE observations, therefore strengthening 
the conclusions above in favour of a core emission origin for the VHE emission of these sources. The model also naturally 
explains the observed very fast variability of the VHE emission. 

Conclusions. 

Key words. Magnetic reconnection - particle acceleration - radiation mechanisms: non-thermal 


1. Introduction 


The non-thermal multi wavelength emission from active 
galactic nuclei (AGNs) has been broadly studied. Regarding 
the very high energy (VHE) emission, until recently only 
AGNs with highly beamed jets towards the line of sight, 
namely blazars, were detected by gamma-ray telescopes. 
More than a chance coincidence, these detections are con¬ 
sistent with the conventional scenario that attributes the 
VHE emission of these sources to particle acceleration along 
the jet being strongly Doppler boosted and producing ap¬ 
parently very high fluxes. 


Lately, however, a few non-blazar sources which belong 
to the branch of low luminosity AGNs (or simply LLAGNs) 
for having bolometric lumino sities of only a f ew times the 
Eddin gton luminosity, LEdd (|Ho et al.||1997t |Nagar et al.| 
have been also detected at TeV energies by ground 
based 7 -ray observatories (e.g., Sol et al. 2013 and refer¬ 
ences therein). The angular resolution and sensitivity of 
these detectors are still so poor that it is hard to establish 
exactly the location of the e mission, i.e. whether it c omes 
from the jet or the core (e.g., Kachelriess et al.||2QlQ ). 


Among these sources, the radio galaxies M87, Centau- 
rus A (Cen A or NGC 5128), Persus A (Per A or NGC 
1275) and IC 310 are probably the most striking cases. 
These VHE detections were surprising because, besides be¬ 
ing highly under luminous, the viewing angle of the jets 


of these sources is of several degrees, therefore allowing 
for only moderate Doppler boosting. These characteristics 
make it difficult explaining the VHE of these sources adopt¬ 
ing the same standard scenario of blazars. 

Eurthermore, observations by MAGIC, HESS and VER- 


ITAS of f 
of IC 310. 

short time scale variability in the 7 -ray emission 

, M87 and Per A (| Aharonian et al.||2006 

|Abdo et 

aLl|2009a| 

i lAckermann et al.||2012| |Aleksic et al. 

[ 2010 a|bl 


very 

In the case of Cen 


it is produced in a 
compact region that might be the core 
A, though there is no evidence of significant vari ability at 
E >lQ0 MeV, GeV or TeV bands by Eerm i-LAT(Ahdo et 
|al.| 2 Q 10 ) or HESS ( Aharonian et al.|20Q9 ), it has been also 
argued that the HESS data of this source wo uld be more 
compatible with a point source near the core (Kachelriess 
|et akpoTO ). If the y — ray photons were due, for instance, to 
proton-proton {pp) interactions along the jet then on leav¬ 
ing the source they would interact with the extragalactic 
background light (EBL) resulting in a flatter spectrum in 
the Te V range which is not com patible with HESS measure¬ 
ments ( Kachelriess et al.|2009b see however, other potential 
explanations in §. 5). 

Though a number of works have attempted to explain 
these observations as produced in th e large scale jets of 
these sources (e.g., Stawarz et al. 2006), the evidences above 
led to the search for alternative particle acceleration sce¬ 
narios involving the production of the VHE in the sur- 


Article number, page 1 of 16 







































A^A proofs: manuscript no. manuscript 


rounds of the nuclear black hole (BH), for instance, in a 
pulsa r-like cascade mechanism in t he BH magnetosphere 
(e.g., Neronov fc Aharonian [2QQ7|), or in the sub-parsec 
ITavecchio & Ghisellinil 


scale j et (e.g 


|2009b). In particular, Tavecchio & G 


20081 

tiisellim 


yoked a two zone model with a jet with a fast spine and 


Abdo et ah 
( |2008p 


m- 


a slow er layer to explain the TeV flares, while Lenain et al. 


(2008) proposed that the emission would occur while the 


jet is collimating, and Georganopoulos et al. (2005) while 
it is decelerating. Another process to explain th ese VHE 
flares has been proposed by Giannios et al. (2010) in which 
misaligned mini-jets driven by magnetic reconnection are 
moving within the jet with relativistic velocities relative to 
it. A two-ste p acceleration model to TeV energies was also 
proposed by 


trying to reconstruct the observed SEDs, specially at the 
VHE branch, of the four radio galaxies mentioned earlier, 
i.e., Gentaurus A, Per A, M87 and IG 310 applying the same 
acceleration model above. 

We first compute the power released by fast magnetic 
reconnection in the surrounds of the BH as described, and 
then the resulting particle spectrum of the accelerated par¬ 
ticles in the magnetic reconnection site. In particular, we 
explore the first-order Eermi acceleration process that may 
occur within the current sh eet as proposed in |de Gouveia 
Dal Pino & Lazarian (2005). Such acceleration mechanism 
has been extensively tested numerically in 3D collisional 
MHD sim ulations of magnetic current sheets employing test 


BH involving initial particle acceleration within the accre¬ 
tion disk and then further centrifugal acceleration in the 
rotating magnetosphere. 

In this work, we consider an alternative model in which 
particles are accelerated, through a first-order Eermi pro¬ 
cess, in the surrounds of the BH by the magnetic power 
extracted from fast magnetic reconnection events occurring 
between the magnetosphere of the BH and the magnetic 
field lines arising from the inner accretion disk (see Eig- 
ure[^. Inspired by similar phenomena occurring in space 
environ ments, like the earth magnetotai l and the solar 
corona, |de Gouveia Dal Pino & Lazarian (2005) explored 
this process hrst in the framework of m icroquasars and then 


particles (Kowal et al 


Istomin fc Sol ( 2009 ) in the surrounds of the al.|2014[) and also in comsloiiless particle in ce l l simulations 


l[2M|[i 

ollisionl 


2012 de Gouveia Dal Pino et 


(e-g-, 

Drake et al.| 2006 |Zenitani et al. 

2009 

Drake et al. 

2010 

Cerutti et al. 2013| 2014| see also the reviews by de 

Gouveia Dal Pino et al.|2014t de Gouveia Dal Pino & Kowal 


2015). 


de Gouveia Dal Pino et al. (2010a) and de Gouveia Dal Pino 
et al. ( |2010bp extended it to AGNs 


In these works, de Gouveia Dal Pino fc Lazarian (2005) 
and de Gouveia Dal Pino et al. (|2010a|) found that fast re- 
connection could be efficient enough to produce the core 
radio outbursts in mic roquasars and AGNs. M ore recently. 


Kadowaki et al. (2015) (henceforth KGS15) and Singh et al. 
(|2015p revisited this model exploring different mechanisms 
of fast magnetic reconnection and accretion, and extend¬ 
ing the study to include also the gamma-ray emission of 
a much larger sample containing more than 230 sources. 
They confirmed the earl ier tr end found by |de Gouveia Dal 


Pino fc Lazarian (|2005|) and de Gouveia Dal Pino et al. 
( |2010ap , an d verified that if fast recon nection is driven by 


turbulence (Lazarian & Vishniac 1999) there is a correlation 
between the fast magnetic reconnection power and the BH 
mass spanning 10^^ orders of magnitude that can explain 
not only the observed radio, but also the VHE luminosity 
from microquasars and LLAGNs (involving all the sources 
of the so ca lled fundamental plane of BH activity [Merlonij 


et al. 


2003). 


The cor relations found in the works above (specially in 
KGS 15 and Singh et al.|2015 ) between the calculated power 


released by magnetic reconnection in the surrounds of the 
BH and the observed radio and gamma-ray luminosities of 
a very large sample of LLAGNs and microquasars, have 
motivated further investigation to test the validity of the 
model. In particular, in a recent work, we explored in de¬ 
tail the non-thermal emission of the microquasars Gygnus 
X-1 and Gygnus X-3 and found that this reconnection ac¬ 
celeration model is able to reproduce most of the features 


In order to reconstruct the SED, we consider the rele¬ 
vant radiative processes due to the interactions of the ac¬ 
celerated particles by magnetic reconnection with the sur¬ 
rounding radiation, matter and magnetic fields. We then 
compare the rates of these radiative losses with the mag¬ 
netic reconnection acceleration rate and determine the max¬ 
imum energy that the electrons and protons can attain. Eor 
comparison, we also consider the acceleration rate due to 
shocks in the surrounds of the reconnection region, but find 
that this is less effective than the acceleration by magnetic 
reconnection. 

We show for the first time that a consistent and numer¬ 
ically tested acceleration model by magnetic reconnection 
in the surrounds of BH sources can effectively reproduce 
the observed SEDs of the four radio galaxies up to TeVs. 

The outline of the paper is as follows. In Section 2, we 
describe in detail our scenario presenting both the accel¬ 
eration model and the emission mechanisms. In section 3, 
we show the results of the application of the acceleration 
and emission model to Gen A, Per A, M87 and IC 310. 
Einally, we discuss and summarize our results drawing our 
conclusions in Section 4. 


2. Description of the acceleration and emission 
model 

As stressed, we consider here that relativistic particles may 
be accelerated in the core of LLAGNs, i.e., in the coronal 
region around the BH near the basis of the jet launching, as 
a result of events of fast magnetic reconnection and examine 
whether this process may reproduce the observed emission 
sp ecially at VHEs. This model has been de scri bed in detail 


by|de Gouveia D al Pino fc Lazarian (2005) andlde Gouveia 
Dal Pino et al. ~ f2Q10ap , and more recently by kGS15. We 


summarize below its main characteristics. We assume that 
the inner region of the accretion disk/corona system alter¬ 
nates between two states which are controlled by changes 
in the global magnetic field. Right before a fast magnetic 
reconnection event, we adopt the simplest possible configu- 
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Coronal 



Fig. 1. Scheme of magnetic reconnection between the lines ris¬ 
ing from the accretion disk into the corona and the lines around 
the BH horizon. Reconnection is made fast by the presence of 
embedded turbulence in the reconnection zone (as indicated in 
the detail). Particle acceleration may occur in the magnetic re- 


connection zone by a first-order Ferm i process (adapted from de 
Gouveia Dal Pino Lazarian||2005 ). 


The magnetosphere around the central BH can be built 
from the drag of magnetic field lines by the ac cretion disk 
(e.g., MacDonald et al.|1986 Koide et al.|2Q02 ). The large- 
scale poloidal magnetic field in the disk corona can in turn 
be formed by the action of a dynamo inside the accretion 
disk or dragged from the surroundings. Under the action of 
disk differential rotation, this poloidal magnetic flux gives 
rise to a wind that partially removes angular momentum 
from the system and increases the accretion rate. This, in 
turn, increases the ram pressure of the accreting material 
that will press the magnetic lines in the inner disk region 
against the lines of the BH magnetosphere thus favouring 
the occurrence of reconnection (see Figure]^. We note that 
according to mean field dynamo theory, an inversion of the 
polarization of the magnetic lines is expected to occur every 
half of the dynamo cycle; when this happens a new flux of 
disk lines should reach the inner region with an inverted 
polarity with respect to the magnetic flux already sitting 
around the BH, therefore, favouring magnetic reconnection 
between the two fluxes. The advection of field lines from the 
outer regions also allows for periodic changes in the polarity. 
Also for simplicity, we assume the co-rotation between the 
inner disk region and the BH magnetosphere. 

The magnetic field intensity in this inner region can be 
determined from the balance between the magnetic pressure 
of the BH magnetosphere and the accretion ram pressure 
and is given by (KGL15): 


B ^ 9.96 X G, (1) 

where m is the BH mass in units of solar mass, ^ is the 
mass accretion disk rate in units of the Eddington rate (^ = 
M/MEdd, with MEdd = 1-45 x 10^®m g s“^), and rx = 
Rx/Rs is the inner radius of the accretion disk in units of 
the BH Schwarzchild radius {Rs = As KGS 15, we 

adopt here rx = 6 . 


2.1. Conditions for fast reconnection in the surrounds of the 
BH 

As discussed in KGS15, the presence of embedded turbu¬ 
lence in the nearly collisional MHD coronal flow of the core 


region of the AGNs can ma ke reconnection very fast (e.g., 
|Lazarian & Vishniac 1999) and cause the release of large 
amounts of magnetic energy power in the scenario described 
in Figure Q 

According t o the model proposed by [Lazarian ik 
[Vishniac^ ( |1999[ ), even weak embedded turbulence causes 
the wandering of the magnetic field lines which allows 
for many independent patches to reconnect simultane¬ 
ously making the global reconnection rate large, Vr ^ 
'^Aiknj^ where Vr is the reconnection 
speed, va is the Alfven speed, and linj and Vturb the in¬ 
jection scale and velocity of the turbulence, respectively. 
This expression indicates that the reconnection rate can be 

as large as Fa. _ 

This theory has been deepl y investigated (e.g. Eyink et 


201 1| Lazarian et al. 


2012 ) and confirmed numerically 


al.. __ 

] Kowal et al.||200^ 2012). In particular, it has been shown 
(Eyink et al.|201l|i that turbulent collisional fast reconnec- 


tion prevails when the thickness of the magnetic disconti¬ 
nuity layer (see Eq. below, and Figure is larger than 
the ion Larmor radius. As demonstrated in KGS 15, for the 
systems we are studying this condition is naturally satis¬ 
fied and we will adopt this process to derive the magnetic 
power released by fast reconnection. We should also notice 
that there has been direct evidences of turbulent reconnec¬ 
tion in solar coronal events 




the Earth magnetotail ( Retino et al.||2Q(T 


Priest 12001) and also in 


The fluids we are investigating here have large hydro- 
dynamical and magnetic Reynolds numbers (KGS 15) im¬ 
plying that they can be easily distorted and become tur¬ 
bulent. For instance, current driven instabilities, can natu¬ 
rally drive turbulence with characteristic velocities around 
the particles thermal speed. Also, the occurrence of contin- 


[ ( 

Liu et al. 

2003) 

w’ 

lich will t 

len be 


niac| (1999) model (see Oishi et al. 2015 [Lazarian et al. 


2015). Numerical simulations of coronal disk accretion also 
indicate the formation of turbulent flow in the surrounds 


of the BH (see e.g, Dexter et al.pOld] ). All these processes 
may ensure the presence of embedded weak turbulence in 
the magnetic discontinuity described in Figure 

The magnetic reconnection power released by turbulent 
driven fast reconnection in the magnetic discontinuity re¬ 
gion (as schemed in Figure]^, has been derived in KGS 15 


^ We note that the strongly magnetized and low dense coronal 
fluid of the systems we are considering in this work satisfies the 
condition L > Imfp > ri (where L is a typical large scale dimen¬ 
sion of the system, Imfp the ion mean free path and ri the ion 
Larmor radius). For such flows a weakly collisional or effectively 
collisional MHD description is more than appropriate and we 
will employ this approach here (see more details in KGS 15). 

^ We note that in the model described here, the turbulence is 
in general sub-Alfvenic due to the strong magnetic fields implied 
and nearly tran-sonic (since the turbulent velocity is of the or¬ 
der of the sound speed and smaller than the Alfven speed) and 
therefore, incompressible. This regime of turbulen ce has been 
extensively investigated in the lite rature (see e.g., [Lazarian fc 


Vishniac|1999| Lazarian et al.|2012 ) and the acceleration formu¬ 
lae employed here have been obtained directly from numerical 
MHD simulations with particle tests injected in current sheets 
with embedded turbulence also implying this regime (Kowal, de 
Gouveia Dal Pino, Lazarian 2012). 
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and is given by: 

tr-1.66x10^ergs“t ( 2 ) 

where I = L/Rs is the height of the corona in units of Rs] 
lx = Lx/Rs^ Lx < L is the extension of the magnetic 
reconnection zone (as shown in Figured see also Table 1), 
q = [1 — and va = is the relativistic 

form of the Alfven velocity, with va^) = B j^ B be¬ 
ing the local magnetic field, p ricrrip the fluid density in 
the corona, ric the coronal number density, rrip the proton 
mass, and ^ = [1 A- in this work, vaq ^ c. The 

results of KGS 15 (Figure 5 in KGS 15) have shown that 
accretion rates ^ between 0.05 < ^ < 1 are able to pro¬ 
duce magnetic reconnection power values which are larger 
than the observed luminosities of LLAGNs. We adopt here 
^ 0.7, but we should notice that the results are not very 

much sensitive to the choice of this parameter. As demon¬ 
strated in KGS 15 and SKG15 studies, one can match the 
observations by taking alternative fiducial combinations of 
the free parameters in the equation above, particularly by 
constraining the size of the height of the corona, L. 

We will employ the equations above in Section 3 to 
model the acceleration of particles in the core of Gen A, Per 
A, M87 and IG 310. The acceleration region in our model 
corresponds to the cylindrical shell around the BH where 
magnetic reconnection takes place (see Figure]^. This shell 
has a length lx, with inner and outer radii given by Rx 
and Rx + ^Rx respectively, where ARx is the width of 
the current sheet given by (KGS 15): 

ARx = cm. (3) 

In §. 3, we will also need the accretion disk temperature 
in order to evaluate the black body radiation field: 


Td = 3.71 X K, (4) 

where 0.05 < a < 1 is the Shakura-Sunyaev disk viscosity 
parameter which we here assume to be of the order of 0.5. 

In addition, we will need the coronal number density to 
compute the radiative losses of accelerated particles which 
is (KGS15) 

Tic ^ 8.02 X cm-^, (5) 


while the coronal temperature is given by 
T, ^ 2.73 X K. 


(6) 


The magnetic power in Eq.[^ heats the surrounding gas 
and accelerates particles. As in KGV15, we assume that 
approximately 50% of the reconnection power is used to 
accelerate the particles. This is consistent with plasma lab¬ 
oratory exp eriments of particle ac celeration in reconnection 
sheets (e.g., Yamada et al. 2014) and also with the o bser- 
vations of flares in the Sun (e.g., Lin fc Hndsnn|p^7T ). 


2.2. Particle acceleration due to magnetic energy released by 
fast reconnection 

A first-order Fermi acceleration may occur when particles 
of the fluid are trapped between the two converging mag¬ 
netic flux tubes moving to each other in the magnetic re¬ 
connection discontinuity with a velocity Vr. |de Gonveia 
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|Dal Pino fc Lazarian (2005) first investigated this process 
analytically and showed that, as the particles bounce back 
and forth undergoing head-on collisions with magnetic fluc¬ 
tuations in the current sheet, their energy increases by 
< AEjE >r-j SVr/3c after each round trip, which leads 
to an exponential energy growth after several round trips. 
If magnetic reconnection is fast, Vr is of the order of the 
local Alfven speed Va and, at the surroundings of relativis¬ 
tic sources Vr va — c and thus the mechanism can be 
rather efficient. 

From the results of 3D MH D numerical simulations of 
this process (Kowal et al.||2012|), we find that the accelera¬ 


tion rate for a proton is given by (KGV15): 


= 1 3 X f T 

^acc,rec,p ^ -lvj y 


- 0.4 


( 7 ) 


where E is the energy of the accelerated proton, Eq = mpC^ ^ 
mp is the proton rest mass, to = lacd'^A is the Alfven time, 
and lace is the length scale of the acceleration region. 
Similarly, for the electrons one can get: 



( 8 ) 


where is the electron rest mass. 

The two equations above do not take into account the 
effects of radiative losses upon the accelerated particles. 
They will be used to compute the acceleration rates in our 
model. 

As stressed in de Gonveia Dal Pino & Lazarian (2005), 
it is also possible that a diffusive shock may develop in 
the surrounds of the magnetic reconnection zone at the jet 
launching region caused by plasmons or coronal mass ejec¬ 
tions. As in solar flares, these can be produced in the recon¬ 
nection layer and released along the magnetic field lines. In 
this case, we expect the shock velocity to be predominantly 
parallel to the magnetic lines, and the acceleration rate for 
a particle of energy E will be approximately given by (e.g., 
Sprmt|p~988 ): 


t 


-1 

acc,shock 


rjecB 
E ’ 


( 9 ) 


where 0 < 77 ^ 1 characterizes the efficiency of the accel¬ 
eration. We fix 77 = 10“^, which is appropriate for shocks 
with velocity Vg ~ 0.1c commonly assumed in the Bohm 
regime. 

In §. 3 we compare both the magnetic reconnection ac¬ 
celeration time (Eqs. [^ & |^ and the shock acceleration 
time (Eq. 1^ with the relevant radiative cooling process 
that cause the loss of energy of the accelerated particles and 
constrain their maximum energy. These particles may lose 
energy via interactions with the surrounding magnetic field 
(producing synchrotron emission), the photon field (pro¬ 
ducing inverse Compton, synchrotron-self-Compton, and 
photo-meson interactions), and with the surrounding mat¬ 
ter (producing pp collisions and relativistic Bremsstrahlung 
radiation). In §. 2.4, we discuss the relevant radiative loss 
processes for electrons and protons which will allow the con¬ 
struction of the SED of the sources M87, Gen A, Per A and 
IC 310 for comparison with the observations. 
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2.3. Particle energy distribution 

The accelerated particles are expected to develop a power 
law spectrum. Their injection and cooling will occur mainly 
in the coronal region around the black hole (see Figure!^. 
The isotropic injection function (in units of erg“^cm“^s^) 
is given by (see e.g. KGV15): 

Q{E) = QoE~Pexp[-E/Emax] (10) 


with p > 0 and E^ax the cut-off energy which can be calcu¬ 
lated by the balance of acceleration and the energy losses. 
Particles can gain energy up to a certain value E^^ax for 
which the total cooling rate equals the acceleration rate. 

We assume for the power law index p values between 1 
and 2.5 for the sources here investigated (see §.3), which 
are compatible with the analytical and numerical studies of 
parti cle accelerat i on within magnet ic reconnection sheets 
(e.g., Drury|[2M2 Kowal et al.||2Q12 ). 

The normalization constant Qo is calculated from the 
total power injected in each type of particle 


L{e,p) = [ (fr f dE E Q(^e,p){E) 

J V j Ern.n.n. 


( 11 ) 


where V is the emission volume in the surrounds of the 
magnetic reconnection acceleration region in Figure (see 
§. 2.5) and p) is the fraction of the magnetic reconnec¬ 
tion power that accelerates the electrons and protons cal¬ 
culated from Eq. This injected power is equally shared 
between protons and electrons. 

The kinetic equation that describes the general evo¬ 
lution of the particle energy distrib ution N(E, t) 


IS 


the Fokker-Pla nck differential equation (Ginzburg & Sy- 
rovatskii 1995). We here use a simplified form of this equa¬ 
tion. We employ the one-zone approximation to find the 
particle distribution, assuming that the acceleration re¬ 
gion is spatially thin enough, so that we can ignore spatial 
derivatives in the transport equation. Physically, this means 
that we are neglecting the contributions to N{E) coming 
from other regions than the magnetic reconnection region in 
the inner accretion disk/corona zone in the surrounds of the 
BH. We consider a steady-state particle distribution which 
can be obtained by setting ^ = 0 in the Fokker-Planck 
differential equation, so that the particle distribution equa¬ 
tion is 


IT^ poo 

NiE) = \ — \-^j^ Q{E)dE. 


( 12 ) 


Here — ^ = Et~^^i , where is the total cooling rate that 
can be calculated assuming all the cooling mechanisms (we 
describe them in the following section briefly). 


2.4. Photon Absorption 

We consider two main absorption processes of the photons 
produced by the accelerated particles in the nuclear region 
of the sources: the gamma-ray photon absorption due to 
e“e+ pair creation, and the absorption of optical and X- 
ray photons due to external interstellar neutral gas and 
dust (photon-neutral) absorption. 


2.4.1. Photon-photon ( 77 ) annihilation 

The produced 7 -rays can be annihilated by the surrounding 
radiation field via electron-positron pair production, i.e., 
7 + 7 ^ e+ + e“. In our model the dominant radiation field 
for this process in the coronal region is due to the scattered 
photons from the accretion disk (see Figure!^. 

To evaluate the optical depth du e to this process, w e 
have adopted the model described in|Gerutti et al.| (j2011|), 
assuming that the 7 -rays are produced within a spherical re- 
gion around the disk with radius extending up to L ^ 20Rs- 
The attenuated 7 —ray lu minosity L^(E^) at a distance z 
above the disk is given by ( [Romero fc Christiansen|[2QQ5| ) 

L^(E^) = (13) 

where is the intrinsic coronal gamma-ray luminosity and 
r{z,Ej) is the optical depth. The calculated optical depth 
depends on the 7 -ray energy and the distance above the 
disk 2 ). 

Figure depicts the gamma-ray absorption spectrum 
for the sources investigated here, for different heights z. 
We see that in all cases at distances larger than ^ O.lRs- 
from the disk surface, the absorption of 7 -rays becomes 
negligible. Since we are adopting here an emission region 
with an extension O.SRs to 20Rs^ it is reasonable to 
exclude the absorption effect above in our calculations of 
the SEDs (see §.3). 


2.4.2. Photon-neutral {'jN) interactions 

The low energy photons produced in the nuclear emis¬ 
sion region will propagate in the surrounding interstellar 
medium of the host galaxy filled mainly by hydrogen and 
helium gas. Photons with energies larger than the hydrogen 
Lyman threshold (13.6 eV) will be able to photo-ionize the 
neutral gas. 

The optical depth resulting from these interactions is 
approximately given by 


Tjh{E^) = Nno-jNiE^) 


(14) 


where Nh is the neutral hydrogen column density, an d 
is the absorption cross section . As in [Reynoso et al.|[ 2 QlTj 
we take this from Ryter ( 1996|) for E^ < 1 keV considering 
that atomic hydrogen and galactic dust are the dominant 
components of the environment. The values of Nh for each 
source investigated here are taken from the observations 
and are listed in Table [H 

This absorption has been considered in the recon¬ 
struction of the SEDs of the four LLAGNs in §. 3. As we 
will see there, the photons produced in the optical-soft X- 
ray range are fully absorbed by these interactions. 


2.5. Radiative cooling processes 

In this section we discuss briefly the relevant radiative loss 
processes for electrons and protons. 

We take into account both leptonic and hadronic ra¬ 
diative loss mechanisms in the emission region. This corre¬ 
sponds to the torus with volume V that encompasses the 

^ We should remark that we found the contribution due to the 
coronal radiation field itself to be much smaller than that of the 
accretion disk. 
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Fig. 2. Left panels: Spectrum of 7 -ray absorption at selected heights z above the plane of the disk in Cen A, Per A, M87, and IC 
310 (a, c, e, and g panels, respectively). Rigt panels: The transmitted flux exp{—T) for different 7 -ray energies as function of the 
height z above the disk in Cen A, Per A, M87, and IC 310 (5, d, /, and h panels, respectively). 
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Table 1. Model parameters for Cen A, Per A, M87 and IC 310. 



Parameters 

Cen A 

Per A 

M 87 

IC 310 

B 

Magnetic field (G) 

1.25 X 10^ 

4812 

1620 

8874 

w 

Magnetic reconnection power (erg/s) 

1.2 X 10^3 

8.2 X 10^3 

5.25 X 10^^ 

4.5 X 10^3 

ARx 

Width of the current sheet (cm) 

3.6 X 1013 

2.4 X lOi^ 

1.35 X 1013 

3.2 X IO 12 

Tie 

Coronal particle number density (cm“^) 

7.1 X 10^ 

10® 

3.3 X 10® 

lOii 

Td 

Temperature of the disk (K) 

1.9 X 10® 

3 X 10® 

5.2 X 10® 

2.25 X 10® 

Rx 

Inner radius of disk (cm) 

8.8 X IQi® 

6 X lOi^ 

5.3 X 1013 

1.7 X lOi^ 

Lx 

Extension of the reconnection region (cm) 

1.5 X IQi^ 

lOi® 

4.4 X 1013 

8.8 X IO 12 

L 

Extension of the corona (cm) 

3 X IQi^ 

2 X 1013 

4.4 X 1013 

8.8 X 1013 

V 

Volume of emission region (cm^) 

7.8 X 10^3 

2.3 X 10^3 

1048 

1.36 X 10^1 

d 

Distance of the source(Mpc) 

3.8 

75 

16.7 

78 

m 

Mass of BH (Mq) 

5 X lOi" 

3.4 X 10® 

3 X 10® 

10® 

p 

Injection spectral index 

2.4 

2.15 

2.4 

1.7 

^min 

Particle minimum Lorentz factor 

6 

2 

4 

2 

Nh* 

* 'T'U. 

Dust/neutral gas column density (cm“^) 

1023 

n 

4 X lO^i 

2 X 10^® 

1.2 X 1031 


Morganti et al. (2008); Canning et al. (2010); Lieu et al. (1996) and Kalberla et al. (2010), respectively. 


cylindrical shell where magnetic reconnection particle ac¬ 
celeration takes place in Figure Considering that the 
cylinder extends up to L in both hemispheres, then the 
small radius of the torus is r = L/2 and the large radius 
is Rx, so that the effective emission zone in our model has 
an approximate volume V = 7t^L‘^Rx- 

For leptonic processes, we consider the interactions 
of relativistic electrons with the surrounding magnetic, 
charged matter, and photon fields. 

Accelerated electrons spiralling in the magnetic field 
emit synchrotron radiation. We calculate the synchrotron 
loss rate for the sources considered here and the radiated 
synchrotron spectrum as functions of the scattered photon 
energy (see Eqs. 10 and 11 in KGV15). Electron interac¬ 
tions with the electrostatic field of nuclei of charge Ze al¬ 
low for the production of bremsstrahlung radiation. Einally, 
relativistic electron interactions with photons may produce 
inverse Compton (IC) radiation. We considered different 
photon fields in the surrounds of BH, namely, the scattered 
photons from accretion disk and core black-body radiation 
(Eq. 36 in KGV15) and the electron synchrotron emission 
(which allows for SSC; Eq. 14 in KGV15) and found that 
the latter, i.e., the SSC mechanism is dominant in the in¬ 
ner coronal/accretion disk region we are interested in this 
work. 

Accelerated protons produce hadronic emission from in¬ 
teractions with the magnetic field (synchrotron), and also 
through the decay of neutral pious. These are produced 
either by inelastic collisions with nuclei of the corona {pp 
interactions; Eqs. 19 & 21 of KGV15), or via interactions 
with photons, in photo-meson production (pq mechanism; 
Eqs. 40 & 43 in KGV15). The latter mechanism takes place 
for photon energies greater than Et^ ~ 145 MeV in the ref¬ 
erence frame of the relativistic proton. A single pion can 
be produced in an interaction near the threshold and then 
decay giving rise to gamma-rays. In our model the domi¬ 
nating photon field comes from the synchrotron radiation. 

E 


^ We find that for photo-meson production, the radiation from 
the accretion disk and from the corona are irrelevant compared 
to the contribution from the synchrotron emission. 


The cooling rates due to these leptonic and hadronic 
mechanisms are plotted as functions of the particle energy 
for the sources here studied in Eigures and [9a| 

and Eigures [31^ [51^ [71^ and [91^ respectively, where they are 
also compared with the acceleration rates due to shock and 
magnetic reconnection. 

The radiative losses considered above are used to calcu¬ 
late the SEDs of the sources in §3. 


3. application to radio galaxies 

We describe here the results of the application of the model 
described above to Cen A, Per A, M87 and IC 310 which 
are classified as radio galaxies. As remarked in §. 1, these 
radio galaxies have been observed at VHE by FERMI- 
LAT, VERITAS and HESS (e. g.,|Abdo et al.|20Q9^[M0 


Abramowski et aL]|2Q12 lAleksic et al.||2Q14bp. 

In general lines, in all cases, the procedure to calculate 
the SED begins with the determination of the total power 
released by fast magnetic reconnection within the accelera¬ 
tion region (Eq. [^, that is, in the cylindrical shell of height 
Lx and thickness ARx (Eq. (see Eigure 1). This power 
is then employed to compute the spectrum of accelerated 
electrons and protons (Eq. that will be injected into 
the emission volume E, i.e., the torus that surrounds the 
acceleration region (as described in §. 2.5). The parame¬ 
ters employed for each source are given in Table 1. We note 
that our model has actually only 7 free parameters (i.e., 
Rx, Tx, L < p^ ^mini ^ and a). The remaining quan¬ 
tities of Table 1 are obtained directly from these parame¬ 
ters through Eqs. 1 to 5 (i.e, 5, IF, ARx, ric and Td), or 
from the observations (i.e, d and m). The maximum energy 
that each particle spectrum can attain is obtained from the 
comparison of the acceleration rates, for both electrons and 
protons, with the relevant radiative loss processes (seen §. 
2). As remarked, the accelerated electrons will loose en¬ 
ergy by synchrotron, IC and Bremsstrahlung mechanisms, 
with a dominance of the synchrotron process shaping their 
spectrum. The fluxes of these emission processes are then 
calculated and also the number density of the synchrotron 
photons that are partially self-scattered by the electrons 
(leading to SSC emission) and by protons (in pq interac¬ 
tions). Likewise, the energy distribution of the protons is 
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also calculated taking into account the radiative cooling 
mechanisms due to synchrotron, pp and pj interactions (see 
§2.3) that will shape the very high energy part of the SED. 


3.1. Application to Cen A 


The Prominent radio galaxy Cen A (or NGC 5128) i s the 
nearest FR I active radio galaxy to Earth (z=Q.Q018, |Gra- 


ham 


1978), at a distance of 3.8 Mpc (Rejkuba |2QQ4p, 


making it uniquely observable among this class of objects 
and an excellent source for studying the physics of relativis¬ 
tic outflows as well as of the core region. Gen A is one of 
the best well known extragalactic objects over a wide range 
of frequencies and the photon emission from the nuclear re¬ 
gion of the galaxy has been detected from the radio to the 
7 -rays band. Gen A has been proposed as a possi ble source 
of UHE cos mic rays (with energies < 6 x lO^^eV; [Abraham 
et al.|[2007 ) by the Pierre Auger collaboration. The SMBH 
mass inferred from kinematics of stars, as well as H 2 and 
ionized gas is estimated to be in the range of r\j 10 ^- 10 ® Mq 

and here we 
angle of the 


2006 


Neumayer et al. 2007) 

lO^M 


(jMarconi et al. 

adopted the value 5 x IO^Mq. The viewing 
jet (0) is still debatable, for instance at parsec scales it is 
6» ~ 50° - 80° (iTingay et al.lHO^, whereas at the 100 pc 
scale 0 ^ 15° (Hardcastle et al. 

In this section we show the results for Gen A obtained 
by applying the model described in §. 2 around the nuclear 
region, employing the set of parameters listed in Table 1. 

The values for the first five parameters in Table 1 have 
been calculated from Eqs. m and We take for the ac¬ 
cretion disk inner radius the value Rx = for the ex¬ 

tension Lx of the reconnection region (see Figure [ 3 , we 
consider the value Lx — lORs, and for the extension of the 
corona L 2 OR 5 . As remarked earlier, the volume V of 
the emission region in Table 1 was calculated by consider¬ 
ing the torus that encompasses the reconnection region in 
Figure The magnetic reconnection power W is evaluated 
from Eq. 

Figure shows the radiative cooling rates for the dif¬ 
ferent energy loss processes for electrons and protons as 
described in Section 2.4. These are compared with the ac¬ 
celeration rates due to first-order Fermi acceleration both 
within the magnetic reconnection site (Eqs.[^&[^ and be¬ 
hind a shock (Eq. 1^. We notice that at hi^i energies for 
both protons and dnctrons the acceleration is dominated 
by the first-order Fermi magnetic reconnection process in 
the core region. Besides, the main radiative cooli ng p rocess 
for the electrons is synchrotron radiation (Figure 3a), while 
for protons the photo-meson production (py interactions) 
governs the loss mechanisms (Figure |3b| ). For the pj in¬ 
teractions, we have found that the proper target radiation 
field is that of the photons from the electron synchrotron 
emission. 

The intercept between the magnetic reconnection accel¬ 
eration rate and the synchrotron rate in Figure [3^ gives the 
maximum energy that the electrons can attain in this ac¬ 
celeration process, which is ^ 3 x lO^^eV. Protons on the 
other hand, do not cool as efficiently as the electrons and 
can attain energies as high as - 2.5 X lO^^eV. 

We have constructed the SED of Gen A using a lepton- 
hadronic model where particles are accelerated close to the 
central BH by magnetic reconnection and interact with the 
surrounding fields radiating in a spherical region of radius 


Article number, page 8 of 16 




Fig. 3. Acceleration and cooling rates for electrons (a) and pro¬ 
tons (b) in the nuclear region of Cen A. 


L (see section 2.4). The SED is depicted in (Figure]^ where 
it is compared with the observations. 

We note that the data presented in this figure in the low 
and intermediate energy ranges come mostly from archival 
data and represent typical average emission and activity. 
The emission in the high energy range obtained by Fermi- 
LAT and HESS at the same epoch correspond to non¬ 
variable and moderate activity too. 

As described in section 2.3, we have adopted a particle 
energy injection power law function: 

Q{E) oc E-P. (15) 

In Figure]^ we considered injected primary particles with 
p = 2A which is consistent with theoretical predictions of 
particle acceleration within magnetic reconnection sites (see 
§. 2.3). Our calculations show that synchrotron radiation 
explains the observed emission in the radio to visible band, 
while SSG is the dominant mechanism to produce the ob¬ 
served hard X—rays and low energy 7 -rays as a result of 
interactions between energetic electrons with scattered syn¬ 
chrotron photons. Also in the Figure, neutral pion (tt^) de¬ 
cays can explain the observed 7 -rays at TeV energies, via 
pp and jry interactions which are the two main processes 
producing tt^. 
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Suzaku ([Markowitz et al.|2006 ). We al so include data trorn OSSE 


( Kinzer et al.|1995| Sind COMPTEL ( Steinle et al.|1998 ) in the 
range of 5 x 10^ — lO^eV. The d ata observed in the energies 


10^ —lO^^eV are taken by EGRET (Sreekumar et al.|l999 


man et aLll999[ 

(Abdo et al.||2QQ9b| [MO) 


Hart- 


(Aharonian et ah |2QQ9p. 


)|) and in the energies' 10° — 10^°eV by' Eermi-LAT 
009bl |2010P. The TeV data are taken by HESS 


We note that in order to fit the observed data in the 
radio to optical range, we had to assume a minimum en¬ 
ergy for the injected electrons in the acceleration zone (Eq. 
11 in KGV15), Emin = {'Imin “ l)TOeC^, with -ymin = 6 . 
Though this injected value has no influence on the VHE 
tail of the SED, it is determinant in the match of the low 
energy branch. We have found that values of jrnin < 6 do 
not lead to the synchrotron match in the low energy range 
(see also §4). 

As remarked in §. 2.4, the 7 -ray absorption due to pair 
production occurs according to Eigures and very near 
the accretion disk at heights smaller than ^ O.OOli^^, thus 
much smaller than the emission region that extends up to 
~ 20Rs in our model, so that exp{—T) 1 and the absorp¬ 
tion effect is not effective at the heights of interest. On the 
other hand, due to the high dust and neutra l gas column 
de nsity i n Cen A {Nh = 10 ^^cm“^, see e.g., Morganti et 


al. 


2008) we find that the optical to soft X-r^ emission is 


y absorbed via jN absorption (see figure 0 Eq. 14) 


3.2. Application to Per A 

Perseus A (also known as NGC 1275 and 3C 84), is a nearby 
active galaxy located at the centre of the Perseu s cluster 
and hosts a central SMBH mass of ^ 3.4 x IO^Mq (Wilman 
et al.||1994). In fact. Per A is one of the closest 7 -ray emit¬ 
ting AgNs. I ts distance to the Earth is 75 Mpc (Brown & 


Adams 1 2011 ) and is also of great interest, specially due to its 
proximity, also providing an excellent opportunity to study 
the physics of relativistic outflows. Per A also seems to ex- 


hibit jet precession with an orientation angle ^ 30° — 55 ° 
( Walker et al.|199^ see also Ealceta-Gongalves et al.|2010 ), 
which may be an indication that Per A is the result of a 


merger between two galaxies (jLiu and ChenJ[2007|). It is a 


very bright radio galaxy showing an extendm jet with ER I 


morp hology (e.g., Vermenlen et al.|[l994| 


] \i9M 

kpc (y 


, jButtiglione et al. 
Pedlar et al.|1990| )' 


Fig. 4. Calculated spectral energy distribution (SED) for Cen A 
using the magnetic reconnection acceleration model in a lepton- 
hadronic scenario compared with observations. The data de¬ 
picted in the radio to o ptical energy range (10 ~^eV — leV) is 
from SCUBA at 800 pm dHawarden et al |1993| ), ISO & SCUBA 
at 450 pm and 850 pm ( [Mirabel et al.||1999[); the data in the 
hard x-rays range is from Swift-BAT (jAjello et al. 2009) and 


2010 ) with asy mmetric jets at bo th ^ 
and pc scales ( |Asada et al.||2006 ). 

The parameters of our model for producing the SED of 
Per A are tabulated in Table 1. The first five parameters are 
calculated from Eqs. Si an dm As for Cen A, we have also 
used for the accretion disk inner radius the value Rx = QRs 
and for the extension Lx of the reconnection region the 
values Lx = lOi^s- and L 20i?s'. 

The radiative loss and acceleration rates for electrons 
and protons are compared in Eigurej^ As in Cen A, mag¬ 
netic reconnection is the dominant acceleration mechanism 
over shock acceleration at the high energy branch for both 
electrons and protons and determines the maximum energy 
that the particles can achieve before losing part of it ra- 
diatively. Electrons may be accelerated up to 3 x lO^^eV 
and the main process to cool them is synchrotron. While 
the maximum energy the protons can achieve is lO^^eV and 
photo-meson production (py) is the dominant mechanism 
to cool them. Similarly to Cen A, the dominant photon field 
interacting with the accelerated protons is the synchrotron 
radiation. 

We have constructed the SED for this source employing 
a leptonic scenario (Eigurej^. In this case, the primary par¬ 
ticles were injected with a power law spectral index p = 2.15 
(Eq. 10). The radio spectrum is matched by electron syn¬ 
chrotron emission, with particles injected into the acceler¬ 
ation zone with rest mass energy (i.e., with 7 min = 2). The 
observed X-ray and 7 -ray emission is nearly reproduced by 
SSC occurring in the nuclear region in a spherical region of 
radius L ~ 20Rs, as described in §. 2.4. However, it should 
be noticed that there is a high intensity BATSE data point 
in the ^ 10^ eV that is not matched by the model. Given 
the fact that this source is highly variable, this feature prob¬ 
ably corresponds to a more active state superposed to the 
less active one (represented by the less intense data points 
in the same energy range). 

The observations also indicate that there is a high en¬ 
ergy cut-off around ^ 3 x lO^^eV in this source. In our 
scenario this is due to leptonic emission produced by inter¬ 
actions of high energy electrons with the radiation field pro¬ 
duced by themselves and this cut-off is compatible with the 
maximum energy calculated from the comparison of the re¬ 
connection acceleration rate with the synchrotron loss rate 
in Eig. 

As stressed in §. 2.5, the optical depth for the produced 
7 -rays was also calculated in this case and is shown if Eig- 
uresand 2 d). We note that the 100 GeV 7 -rays may be 
fully absorbed due to pair production but only very near 
to disk {z < O.OOlRs). However, these vertical distances 
from the disk, comparing to the length scale of the emis¬ 
sion region is very small and reasonably, we can ignore the 
absorption effect at the heights larger than ^ Rg which is 
compatible with extension of emission region in our model. 

Similarly to Gen A, the neutral gas and dust of the 
interstellar medium of the host galaxy in this source also 
causes the extinction of the em ission in the range of 10 — 
lO^eV (Eig. see also §. 2.4.2). 
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(a) 



Fig. 5. Acceleration and cooling rates for electrons (a) and for 
protons (b) in the nuclear region of Per A (NGC 1275). 


3.3. Application to M87 


The FR I giant radiogalaxy M87 is another well-known 
nearby AGN located at 16.7 Mpc within the Virgo cluster 
which harbours a SMBH with a m ass of Mbh - 6 x 10 ^ M© 
(e. g., Gebhardt fc Thonias|2QQ9 ) which, along with Cen A 
and Per A, has been known as a peculiar extragalactic lab¬ 
oratory to study high energy astrophysics and investigate 
the nonthermal mechanisms of VHE emission in AGNs. The 
observations ind icate that its jet is oriented within 20 ° of 
the line of sight ( Biretta et al.|1999 ), so that as in the other 
cases, no significant Doppler boosting is expected for the 7 - 
ray flux. 

The T eV 7 -ray signal from M 87 was first reported by 
HEGR A dAharonian et ^J20Q3 ) and then confirmed by 


HESS (Aharonian et al. 


. The latter also revealed that 


this emission is strongly variable with time scales of 1-2 days 
and thus produced in a very compact region, as pointed out 
before. 


Table 1 shows the parameters that we used to calculate 
the acceleration and cooling rates and also to reconstruct 
the SED of this source. 

In Eiguresj^ we compare the rates of the radiative cool¬ 
ing processes with the rates of the acceleration mechanisms 
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Fig. 6. A leptonic model to reproduce the SED of Per A (NGG 
1275) using the ma gnetic reconnectio n accelera tion modeh Dat a 
includ e MOJAVE jhister et al.|20'09 ), Planck ( Ade et al |2011|), 
HST JCtiiaberge et al.| | 1999p , and HST FOS ( [Johnstone ^ 
Fabianll 19951) for the radio t o optical spec trum; data d epicted 
in X- rays is from the XM M (|TQrresi||2012 ), Swift-BAT (Ajello 
et al.|2009), and BATSE (Harmon et al.| 2004); and data de- 
pictea uTThe gamma-ray band is from Fermi-LAT (Abdo et al. 
2009a Ackermann et an|2012|) in the 100 MeV- lOOGeV energy 


range, and from MAGIC (lAleksic et al.|2010a|b| 2012a|b ) in the 
VHE tail. We note that the error bars f or the HATSE data (in 
the 10^ eV range) wer e evaluated u s ing Harmon et al. (2004); 
Soldi et al. (2014) and Wilson et al. (2 (jl2[ ). 


for electrons and protons. As in the other two sources, we 
find that the dominant energy loss mechanisms are the syn¬ 
chrotron and the p 7 interactions for electrons and protons, 
respectively, and the acceleration is dominated by the mag¬ 
netic reconnection process which defines the energy cut off 
for both electrons and protons. Eigure indicates that 
this maximum energy is ^ 4 x lO^^eV for electrons and 
^ 5 X lO^^eV for protons. 

Eigureshows the calculated SED for M87 compared to 
the observations. It is also reproduced by a lepton-hadronic 
model in the core region as described in §. 2 , where we 
assumed an injected particle energy distribution oc E~^ 
with a power index p = 2.4. 

With an electron minimum energy Emin = 4meC^, we 
can fit the observed core radio to visible spectrum by syn¬ 
chrotron emission. 

As in Cen A, the low and intermediate energy data in 
Eigure come from archive and represent typical average 
emission. The data obtained by Eermi-LAT (10^ — lO^^eV) 
and by HESS (the TeV tail) correspond both to more 
quiescent states and have been taken in different epochs. 
They are reproduced in our model by different mechanisms. 
While Eermi-LAT data are fitted by SSC and pp collisions, 
HESS data are fitted by the decay of neutral pious from 
pj interactions with photons coming from the synchrotron 
radiation. We note that an observed flare state by HESS 
(Aharonian et al. 2006; not shown in Eig.[^ can be also re¬ 
produced by our model assuming a flatter injection particle 
spectrum with a spectral index = 2 . 1 ). 

Eigures and show the absorbed 7 —ray flux for 
M87. As in the other cases, this absorption is significant 
only at heights smaller than Rs and therefore, its effect can 
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(a) 



Fig. 7. Acceleration and cooling rates for electrons (a) and for 
protons (b) in the nuclear region of M87. 


be neglected at the much larger emission scales considered 
here. 

The absorption of low energy photons by interstel¬ 
lar neutral gas and dust in this source is also important 

(Fig.|§. 

We note that our model and the chosen parametrization 
is also consistent with the observed TeV rapid variab ility of 
M87 which is ^ 1 — 2 days ( Abramowski et al.|2Q12 ) imply¬ 
ing an extremely compact emission region (corresponding to 
scales of only a few Rs)^ As remarked in §2.4, the emission 
region in our model corresponds to the torus region that 
encompasses the cylindrical shell where magnetic reconnec¬ 
tion particle acceleration takes place in Figure i.e., the 
effective emission zone in our model has a thickness L. 
For this source L bRs (see Table 1), which is of the order 
of the inferred scale from the observed variability. 


3.4. Application to 1C 310 

The peculiar galaxy IC 310 (also named B0313+411 and 
J0316+4119 in observational reports) is one of the brightest 
objects which, as Per A, is also located in the Perseus galaxy 


cluster at a distance of 78 Mpc from Earth (Aleksic et al. 
2014c) and harbours a supermassive BH with a mass 





Fig. 8. A lepton-hadronic model of the SED of M 87 compared 
with observati ons. The core radio dat a are obtained fro m MO- 
JAVE VLBA ( Kellermann et al.|2004 ) at 15 GHz, from (|Biretta 
1^911) at 1.5, 5 and 15 GHz, from IRAM (|Desmingreet 
at 89 GHz, from SMA at 230 GHz ([Tan et al.||2008| ) 



1996). 




7 -ray . 


ronian et al.|2006] ) from HESS. 


^ 10 ^ Mq ( Aleksic et al.| 20 f^ ). The redshift of this source 
is z=0.0189 (|Bernardi et al.||2002p which has made it the 
fourt h nearest AGN at VHE gamma-rays (Kadler et al. 
2012), after Cen A with z=0.00183, M 87 with z=0.004 
and Per A with z=0.017559. 

IC 31 0 has been observed at energies E >100 GeV by 
MAGIC (jMariotti et al. 2010) and Eermi-LAT collabora- 
tion also reported th e detection of photons above 30 GeV 
(Neronov et al.|2010). However the origin of the gamma-ray 
emission is not clear yet and both the jet and the core have 
been considered as possible emission regions. 

Recently, MAGIC collaboration has reported fast time 
variability for IC 310 on the VH E 7 -ray with time scales 
^4.8 min (Aleksic et al. 2014a|c ) which constrains the size 
of the emission zone to 20 % of its Rs- 

The parameters we used to calculate the acceleration 
and cooling time scales and also to reconstruct the SED of 
this source are shown in Table 1. 

The comparison between the acceleration and cooling 
rates are depicted in Eigures 9a and 9b for electrons and 
protons, respectively. As in the other cases, we see that the 
calculated maximum energy for both electrons and protons 
reaches larger values for magnetic reconnection than for 
shock acceleration, so that magnetic reconnection should 
be the dominating mechanism to accelerate particles in the 
nuclear region of this source as well. The diagrams indicate 
that electrons can accelerate up to 8 x 10^^ eV, while the 
protons up to 2 X 10^^ eV. Also in this source synchrotron 
emission is the dominant loss mechanism for electrons and 
P 7 radiation is the dominant one for protons for energies 
larger than r\j 10^^ eV. 

EigurepT] shows the calculated SED for IC 310. As in the 
other sources, the observed radio emission can be explained 
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Log(EQ/eV) 


(a) 



Fig. 9. Acceleration and cooling rates for electrons (a) and for 
protons (b) in the nuclear region of IC 310. 


by synchrotron and the TeV 7 -rays by the pp and py pro¬ 
cesses due to particles injected with a power law spectral 
index p=1.7. The core opacity to this emission has been also 
calculated for IC 310 in figure which indicates that the 
7 -ray absorption is negligible in the emission length scales 
here considered which are above O.SRs- 

Also in this case the radiation in low energy range (10 — 
lO^eV) is fully absorbed due to photon-neutral interactions 
(Figure 

As for M87, our model and the adopted parametrization 
can also naturally explain the fast variability of the VHE 
7 -rays in 3C 310. The effective emission zone for this source 
is L O.SRs according to our model (see Table 1 ), which is 
compatible with the scale inferred from the observed high 
variability in the 7 -ray emission, ^4.8 min. 

4. Discussion and Conclusions 

Our main purpose here was to explore an alternative mech¬ 
anism to explain the spectral energy distribution (SED) 
and, particularly, the VHE emission in non-blazar AGNs, 
i.e., those for which the jet probably does not point to the 
line of sight and therefore, do not have their emission en¬ 
hanced by Doppler boosting in the relativistic jet. Instead, 
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Fig. 10. A lepton-hadronic model of the SED of IC 310 com- 
pared with observations. The core radio data are obtained from 


Kadler et al.||2012| |Dunn et al. 112010 

Becker et al.||19911 |White| 

V Becker||1992[ |Condon et al.|l2002 

Douglas et al.||19961 The 

X-ray data are from XMM-Newton 

Sato et al.|p005|) and the 

VHE 7 -ray data are from MAGIC (I 

deksic et al.|2014c). In the 


upper right side of the diagram it is depicted the detail of the 
modeling of the VHE branch. 


we have examined an acceleration mechanism occurring in 
th e innermost region of the AG N. Based on recent results 


by de Gouveia Dal Pino et al. (2010a), KGS15, and Singh 
et al. ~ pQ15p we investigated the role of fast magnetic re¬ 
connection events in accelerating particles in the nuclear 
regions of low luminosity AGNs, applying this acceleration 
model to reconstruct the SED of Gen A, Per A, M87 and 
IG 310. 

According to this model, trapped particles within the 
magnetic reconnection discontinuity formed by the en¬ 
counter of the magnetic field lines arising from the accre¬ 
tion disk with those of the BH magnetosphere (Eigure 1), 
can be accelerated to relativistic velocities by a first-order 
Eermi process in the surrounds of the B H, as described in 
de Gouveia Dal Pino & Lazarian (2005). Magnetic recon- 
nection events will occur specially when there is substantial 
increase in the disk accretion rate which helps to press the 
two magnetic fluxes together. Since turbulence is expected 
to be present in these systems (s ee §. 2 . 1 ), the reconnection 
can be made naturally fast by it ( [Lazarian fc Vishniac|1999| 
and Kowal et al.||2009 ) 


KGS15 and Singh et al. (2015) used this model to com¬ 
pute the magnetic power released by reconnection in the 
surrounds of the BHs and compared it with the core radio 
and 7 —ray luminosities of outbursts of a sample with more 
than 230 sources including the microqnasar s and LLAGNs 
of th e fundamental plane of BH activity (Merloni et al. 
2003) and also hundreds of blazars and GRBs. They have 
found that the reconnection power is large enough and 
is correlated with the observed luminosities of the micro¬ 
quasars and LLAGNs, following the observed trend for more 
than 10 ^^ orders of magnitude in luminosity and mass of 
these sources. E] 


^ We note that this power is generally not enough to reproduce 
the radio or gamma-ray luminosities of the blazars and GRBs, 
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These correlations found with the emission of micro¬ 
quasars and LLAGNs are very important, because they con¬ 
nect the non-thermal emission, specially the VHE one, with 
the core of these sources, offering a reliable, self-consistent 
mechanism for their origin. However, in order to determine 
the real effectiveness of this mechanism, it is necessary to 
reproduce the SEDs of these sources based on this model. 
This has been done recently for two microquasars Cyg XI 
and Cyg X3, for which the observed SEDs have been suc¬ 
cessfully reproduced by the model above (KGV15). In this 
work we have extended this analysis applying the reconnec¬ 
tion model to the four radio galaxies with detection of VHE 
emission up to TeV. 

Considering all relevant leptonic and hadronic radiative 
loss times due to the interactions of the accelerated particles 
with matter, radiation and magnetic field in the core regions 
of Gen A, Per A, M87 and IC 310, we compared these times 
with the acceleration times and found larger energy cut¬ 
offs for particles being accelerated by mametic reconnection 
than by diffusive shock (see Eigures and0. This 

result stresses the importance of magnetic reconnection as a 
potential acceleration mechanism in the core regions around 
BHs and compact sources in general. 

Moreover, we note that the maximum energies for the 
electrons and protons obtained from these comparisons, are 
actually much smaller than the maximum possible energy 
that the particles can attain within the acceleration region 
in the reconnection sheet. The latter is constrained by the 
thickness of the acceleration region, i.e., ARx (given by 
Eq. which must be larger than or equal to the par¬ 
ticle Larmor radius, vl = EjceB. This implies that the 
maximum energy to which the particles can be accelerated 
by magnetic reconnection is 10^0 eV, SO that the recon¬ 
nection layer is large enough to accelerate the particles to 
UHEs. This value reassures the efficiency of this accelera¬ 
tion process and suggests that if the surrounds of the BHs 
in AGNs were not so full of interacting matter and radia¬ 
tion fields they might be excellent sites for the production 
of UHECRs. 

The cut-off values above were employed in the determi¬ 
nation of the accelerated particle fluxes in the construction 
of the SED of the sources. The electron synchrotron cooling 
rate is found to be larger than any other loss mechanism 
in all leptonic energy range (Eigures 5a, ^ and [9a|). 
Also, it is the dominant process proving tne radiation 
field that produces the SSC and the photo-meson (pq) ra¬ 
diation in the SEDs of the four sources. We have also found 
that the pq process is the dominant mechanism to cool the 
accelerated relativistic protons in the high energy branch. 
This is shown in Eigure [31^ for Gen A for proton energies 
^ 6 X lO^^eV, and in Eigure pm for Per A for proton en¬ 
ergies ^ 2 X lO^^eV, while in figures ^ and 9b for M87 
and IG 310, we see that synchrotron radiation can cool the 
protons for energies ^2xl0^^eV, but for larger energies it 
is overcome by the pq mechanism. 

In summary, we have shown that, employing fiducial 
parameters, our acceleration model is capable of explaining 
the non-thermal low and high energy emission of the SED 
of the four investigated LLAGNs Gen A (Eigu re B , Per A 
(Eigure [^, M87 (Emgre|^ and IG 310 (Eigure [T^. 


but this is compatible with the notion that the emission in these 
cases is produced outside the core, at the jet that points to the 
line of sight and screens any deep nuclear emission. 


An interesting advantage of the model presented here 
is the relatively small number of free parameters used to 
construct the SED (seven). In particular, we adopted an 
one-zone approximation in order to avoid the introduction 
of more free-parameters. Eurthermore, this work is a first 
attempt to test the acceleration model by magnetic recon¬ 
nection in the surrounds of these BH sources whose physics 
and distribution of the photon, density and magnetic fields 
are still poorly known making the use of multi-zone mod¬ 
els e ven more challenging or a rtificial. Unlike other models 
(e.g. Kachelriess et aL||20^09b ), our acceleration mechanism 
was used to constrain the cnaracteristics of the accelera¬ 
tion region. While most of the models take the maximum 
ener gy of the particles as a free parameter to fit the SED 

I |Kachelriess et al.||2010|), 
y from the acceleration 


(e.g. Abdo et al. | |2009a|b|c| |2Q10| 
our model determines this direct 
model, as described. 


4.1. Cen A 


According to our results for Cen A, the observed hard X- 
rays and Femi — LAT q-ray data can be interpreted as 
due to SSC and pp interactions, respectively, with acceler¬ 
ated particles injected in the nuclear region (at distances 
^ 20Rs) driven by magnetic reconnection with a distri¬ 
bution with power law index p = 2.4. The TeV radiation 
observed by HESS, on the other hand, is explained by 
neutral pion decays resulting from pq interactions. 


In Sahn et al.|2012[ the authors also showed that the TeV 
q-rays in Cen A could be explicated by pq interactions, but 
between relativistic protons accelerated by Eermi process 
in shocks along the jet with the mo nochromatic photons 


observed at 170keV. Another model ( Reynoso et al.||2Qll 
also proposed particle acceleration at the jet basis witn the 
production of the hard X-rays by synchrotron emission, and 
the Fermi — LAT and HESS da ta by IC and pp mecha- 
nisms, respectively, along the jet. [Kachelriess et ar]|201Q[ 
on the other hand, have argued against the production of 
the q — rays in Cen A by pp interactions along the jet be¬ 
cause on leaving the source they would interact with the 
EBL resulting in a flatter spectrum in the TeV range than 
the observed one by HESS (see also §. 1). 

All these studies demonstrate that the origin of the VHE 
emission in this source is still highly debatable. A core ori¬ 
gin, as the one suggested in this work arises as an interesting 
possibility as long as magnetic activity is significant in the 
surrounds of BHs, and it might be considered as well. To 
disentangle this puzzle we will need substantially improved 
observations, specially in the q—ray range. It is possible 
that with the much larger resolution and sensitivity of the 
forthcoming CTA observatory (Actis et al. 2011;Acharya et 
al. 2013; Sol et al. 2013), and with longer times of exposure 
of this nearby source (and also of M87, Per A and IC 310), 
we may collect higher resolution data, and more significant 
information on variabiiity that may heip to determine the 
iocation of the emission region. 


4.2. Per A 

In the case of Per A, there is no reievant data yet in TeV 
energies, but our core modei can neariy expiain the observed 
Fermi — LAT and MAGIC data in the 0.1 GeV — 650 
GeV range with a ieptonic scenario dominated by SSC. The 
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synchrotron photons that are absorbed in SSC are produced 
by accelerated electrons by magnetic reconnection in the 
coronal nuclear region around the BH (within distances ^ 
20Rs) having a distribution with a power law index p = 
2.15. 

An SSC model has been also proposed by |Aleksic et al.| 
(2014b), but they assumed that the Per A core could be 
a BL Lac blazar with the jet bending strongly at larger 
scales and the high energy non-thermal radiation could be 
originated in a sub-structure of the jet near the core point¬ 
ing towards our line of sight. This bending still requires 
observational support and any evidence of jet prece ssion 


Pino fc Lazarian||20Q5l de Gonveia Dal Pino et al.||2QlQap . 


However, Giannios et al. (|2010p study provides no predic- 

tions for the SED structure of M87. 

In addition, there is an extensive list of models that 
propose that the variable VHE emission of M87 can be 
originated in the inner jet. These s pan from leptonic mod¬ 


els, s uch as the decele rating flow dGeorganoponlos et al. 


2005), the spine-shear (Ta vecchio fc Ghisellin i 2008), and 
the mini/multi-blob model ( [Lenain et al.|2008( ), to hadronic 
models wit h the emission due to proton synchrotron-py in¬ 
teractions ( Reimer et al.|200'4 ), or pp interactions in a cloud- 
jet scenario ( [Barkov et al.||20f2 ). However, the location of 
the emission region is still an unsolved problem. 


Neronov fc Aharonian (2007) also proposed a nuclear 
origin for the TeV emission of M87 coming dire ctly from 
the magnetosphere of the black hole (see also [Levinson 


2000). They showed that accelerated electrons in the strong 


i 


(e.g.. Walker et al [Falceta-Gongalves et al. 2010| ) may 
favour this model. But our proposed model dismisses the 
necessity of such a strong bending and besides, is supported 
by t he correlations wit h the observations found in KGS 15 
and Singh et al. (2015). 


4.3. M87 

In the case of M87, we have applied the same magnetic 
reconnection model in the nuclear region around the BH 
(within a region of 20 Rs) considering the injection of the 
accelerated particles with a power law index p = 2.4. This 
has resulted a lepton-hadronic scenario for the SED with 
SSC emission and neutral pion decays from pp collisions 
explaining the Fermi — LAT data. We also found that the 
decay of the neutral pious due py interactions can explain 
the observed data by HESS in the TeV range. 

The suggested sites of TeV emission for this source in 
for mer works range from large scale structures of the kpc 
jet (jStawarz et al. 2005|) to a compact peculiar hot spot 
(th e so-called HST-1 kn ot) at a distance 100 pc along the 
jet ( jStawarz et al.|20Q^ and i nner (sub-parsec) parts of the 
sources. Reynoso et al.f(2011) for instance, considered that 
this emission is produced in the jet, but reconstructed all 
the emission features, which are highly variable and possi¬ 
bly non-simultaneous, with a single pp mechanism. 


Giannios et al. (2010), on the other hand, proposed that 
compact minijets induced by magnetic reconnection moving 
relativistically within the jet in different directions, some of 
which pointing to our line of sight, might explain the short- 
time variable TeV flares observed in M87. This model bears 
several similarities with ours as it proposes that the minijets 
are generated by reconnection events in the core region, and 
then move out with the jet flow up to scales of ^ lOOR^. 
Our model also predicts the development of outbursts with 
the formation of reconnected features (plasmons) that may 
be carried out with the jet and might expl ain e.g., observed 
snperluminal features near the jet basis (jde Gonveia Dal 


rotation-induced electric fields in vacuum gaps in the BH 
magnetosphere, similar to a pulsar magnetosphere, could 
lead to the observed TeV emission. Since the acceleration 
and emission mechanisms occurs in a very compact region 
close to the event horizon of the BH, it potentially can ex¬ 
plain the observed variability of TeV y-ray emission from 
M87. Besides, as in our model, they also explain this emis¬ 
sion as due to py interactions with an IR compact target 
photon field produced by synchrotron emission. However, 
as stressed in §. 3, the attenuation of y-ray emission due 
to electron-positron pair production may be significant in 
distances smaller than or equal to ^ Rs (Figure]^, which 
may affect their results. In our model, the emission scales 
are larger (~ ^Rs) making these attenuation effects negli¬ 
gible. 

4.4. ICS 10 

In the case of IG 310, also a lepton-hadronic model in the 
nuclear region is able to explain the observed SED features 
with protons and electrons accelerated by magnetic recon¬ 
nection and injected in the emission region with a power law 
distribution with index p = 1.7. As remarked, the observed 
radio emission is well fitted by synchrotron and the VHE 
emission detected by MAGIC can be explained by decays 
of neutral pious resulting from pp and py interactions. 

Our model with an appropriate choice of parameters is 
also able to explain naturally the time variability detected 
in the sources here investigated. In particular, the very fast 
variability reported for the IC 3 10 y—ray emission of about 
^ 4.8 min ( Aleksic et al.||2Q14c ) implies an emission region 
scale of ^ 0.37^5'. To explain this variability and compact¬ 
ness of the emission region, [Aleksic et al.||2Q14^ suggested 
that the particles could be accelerated by electric fields in 
the BH magnetosphere, as in pulsar models. Nevertheless, 
as demonstrated, the model described here reproduces the 
observed SED with an emission region with a similar size 
as above. 

In conclusion, in the construction of the SEDs for the 
sources discussed here (Gen A, Per A, M87 and IC 310) 
based on our magnetic reconnection model in the core re¬ 
gion, the observed emission at low energies (radio to opti¬ 
cal) can be explained by synchrotron emission. SSC with 
target photons coming from electron synchrotron emission 
is the dominant (leptonic) mechanism to produce the ob¬ 
served hard X-rays and low energy y—rays, while neutral 
pion decays resulting from pp inelastic collisions is the dom¬ 
inant hadronic process to produce the high energy y—rays, 
and neutral pion decays resulting from photo-meson inter¬ 
action (py) the dominant one to produce the very high en¬ 
ergy (VHE) y—rays. Interestingly, in the case of the mi¬ 
croquasars Cyg XI and Cyg X3, KGV15 have also found 
that the core model could reproduce the full observed SED 
including the low and high energy branches. 


® We should remark that the observed emission at low energies 
(radio to optical) from the core regions in the case of Cen A and 
M87 is fitted by the core model described here only if we assume 
that the minimum electron energy injected in the acceleration 
region is a few times the particle rest mass. If one considers 
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4.5. Conclusions 

— We have presented a reconnection acceleration model in 
the core region around the BH of the low luminous radio 
galaxies Cen A, Per A, M87 and IC 310 and showed that 
it is able to reproduce very well their SEDs, from radio 
to gamma-rays up to TeV energies. 

— This is complementary to a recent study where we have 
performed similar analysis for the galactic black hole bi¬ 
naries, i.e., the microquasars Cyg XI and X3 (Khiali et 
al. 2015). Together, these two works strengthen the con¬ 
clusions of the previous works of Kadowaki et al. (2015) 
and Singh et al. (2015) in favour of a core emission spe¬ 
cially for the observed gamma-ray radiation of micro¬ 
quasars and LLAGNs (which belong to the so called 
fundamental plane of BH activity). 

— Magnetic reconnection acceleration seems to provide a 
better efficiency in regions where magnetic activity is 
dominant in comparison with diffusive shock accelera¬ 
tion as the cores of LLAGNs. Particles can gain energy 
up to a few times ^100 TeV due to magnetic reconnec¬ 
tion acceleration. 

— The observed TeV 7 -ray emission may be originated 
in these cores via neutral pion decays in hadronic pro¬ 
cesses. 

— The fast magnetic reconnection acceleration model oc¬ 
curring in the core of these sources can naturally explain 
the observed short time variability, specially of the high 
energy 7 -ray. 


Finally, we note that it is possible that a neutrino 
spectrum may be also produced in the nuclear region of 
LLAGNs considering the same model here investigated, as 
due to charged pion decays via pp and p 7 interactions. In 
a concomitant work, this possibility has been investigated 
to explain the recently obser ved extragalactic neutrino ffux 


by th e IceCube experiment (Khiali & de Gouveia Dal Pino 


2015). 
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